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Oh 

O Abstract. Analytical investigation of time-dependent accretion in disks is carried out. We consider a time-dependent 
disk in a binary system at outburst which has a fixed tidally-truncated outer radius. The standard Shakura-Sunyaev 
I model of the disk is considered. The vertical structure of the disk is accurately described in two regimes of opacity: 
■ Thomson and free-free. Fully analytical solutions are obtained, characterized by power-law variations of accretion 
I rate with time. The solutions supply asymptotic description of disk evolution in flaring sources in the periods after 
^ • outbursts while the disk is fully ionized. The X-ray flux of multicolor (black-body) a-disk is obtained as varying 
' quasi-exponentially. Application to X-ray novae is briefly discussed concerning the observed faster-than-power decays 
fS| ■ of X-ray light curves. The case of time-dependent advective disk when the exponential variations of accretion rate can 
occur is discussed. 

o . 
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1. Introduction 
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-4-^ I We investigate analytically the problem of time-dependent accretion, closely related to the phenomena of flares widely 
^ I observed in binary systems. In this paper we will concern ourselves with the emission of the flaring source which is 
generated by the accretion disk. We suppose the light curve to be regulated by the accretion rate variations. Such 
sources are typified by the low-mass X-ray binaries and cataclysmic variables. 



' After Weizsacker ( 1948 ) who considered the evolution of a protoplanetary cloud, the analytical investigations of 



?H I non-stationary accretion were carried out by Liist ( 1952 ), Lynden-Bell & Pringle ( |I974[ ), Lyubarskii & S haku ra ( |1987| , 



hereafter LS87) as apphed to accretion disks. A brief overview is presented in the book by Kato et al. ( 199^ ) 



LS87 suggested three stages of evolution of a time-dependent accretion disk (see also § 5.1). Initially a finite torus 
of the increased density is formed around a gravitational centre. Viscosity causes the torus to spread and develop into 
the disk (1st stage). After the disk approaches the centre, the accretion rate reaches the maximum value (2nd stage) 
and begins to descend (3rd stage). During this stage the total angular momentum of the disk is conserved. 

Ogilvie ( 1999| ) presented a time-dependent self-similar analytical solution for a quasi-spherical advection-dominated 



flow with conserved total angular momentum. 

In a binary system, variations of accretion rate can be due to the non-stationary exchange of mass between the 
components of the binary (mass-overflow instability model) or due to the disk instability processes (see Kato et 



al. 199S and references therein). At some time the accretion rate onto the centre begins to augment. We assume 
that the maximum accretion rate through the inner boundary of the disk corresponds to a peak of outburst and the 
accretion rate decreases afterwards. 

In this study we particularly focus on the stage soon after the outburst. In § 2 we outline the general equations of 
time-dependent accretion. The basic equation relates surface density of the disk and viscous stresses in it. Thus the 
specific structure of the disk influences greatly the run of the process. § 3 introduces the investigation of time-dependent 
Keplerian a-disks. The vertical structure of standard Keplerian disk is considered in § 4. In a binary system the third 
stage of LS87 cannot be realized because the accretion disk around a primary would be confined by the gravitational 
influence of a secondary. Such disks do not preserve their angular momentum, transferring it to the orbital motion. 
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We suggest the particular conditions at the outer boundary of the disk which aUow the acquisition of new solutions 
characterized by faster decays than in LS87. The procedure and the analytical solution are presented in § 5. 

We calculate the resulting bolometric light curve taking into account the transition between the opacity regimes 
as accretion rate decreases (§ 6). We note that observed light curves can have different slopes due to unevenness of 
spectral distribution (§7). 

In § 8 we discuss the case of advection-dominated accretion flow (ADAF) in which exponential variations with 
time of accretion rate possibly take place. 

In § 9 we discuss application of our model to X-ray novae. 



2. Basic non-stationary accretion disk equation 



In the approximation of Newtonian potential we assume that the velocity of a free particle orbiting at distance r 
around a gravitating object is 



(1) 



where Wk is the Kepler angular velocity; AI is the mass of the central gravitating object, constant in time; G = 
6.67 X 10~* cm^ g~^ s~^ is the gravitational constant. This is a good approximation to the law of motion for particles 
in the standard sub-Eddington disk. In the advection-dominated accretion flow (ADAF) the particles are substantially 
subjected to the radial gr adien t of pressure and thus have the velocity different from that given by (|l|). Following the 
model by Narayan & Yi ( 1994 ), one can assume that the angular velocity in ADAF is = C2 • 

The height-integrated Euler equation on ip and the continuity equation along the height Z are given by: 



So Vr 

dt 



dr 



1 d 

— — EoWrf 

r or 



(2) 
(3) 



where uj is the angular velocity in the disk; Eo(?', ^) = 2 / pdZ - the surface density of the matter, and Wrip{r,t) — 



2 / Wrtp dZ is the height-integrated viscous shear stresses between adjacent layers. The time-independent angular 



velocity is assumed although there can possibly be certain variations of uj in the non-Keplerian advective disks when 
a time-dependent pressure gradient is involved (see, e.g. Ogilvie 199£). 

It is convenient to introduce the following variables: F = Wripv'^, henceforth 2ttF means the total moment of 
viscous forces acting between the adjacent layers, h^, ~ wr^ - the specific angular momentum of the matter in the 
disk, and h = bJ^.r'^ 

M{r,t) 



From Eq. (|) in view of (|l|) it follows that 



2tt 



'dh 



OF 
'dh 



Substituting 



9So _ 1 {GMf d_ 
dt ~ 2 dh 



in (^ and expressing r in terms of h, we obtain the basic equation of time-dependent accretion: 



dh 



(4) 



(5) 



In the case of the Keplerian disk dh^,/dh ~ 1. The advection-dominated solution by Narayan and Yi ( |l994|) yields 
dh^/dh — C2, where C2 is a dimensionless constant. 



3. Non-linear problem of evolution of the standard Shakura-Sunyaev disk 



The special case when the moment of viscous forces depends linearly on the surface density and has a power law 
dependence on the radius {F oc Sq^') was thoroughly investigated by Lynden-Bell & Pringle (1974). In this particular 
case Eq. (|^) is linear and the solution can be presented as the superposition of particular solutions (Green's functions) 
while the non-linear equatio ns do not allow such solut ions. In the paper by LS87 the necessary relation between Eq 
and F for a-disks (Shakura 1972 ; Shakura & Sunyaev 1973 ) was derived from the vertical structure equations. Then 
Eq. (0) acquires the following non- linear form taking into account that h = h^:. 



dF 
'dt 



D 



pm Q2p 

1F 'dh? 



(6) 
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where D is the dimension constant; m = 2/5, n = 6/5 when the Thomson scattering dominates the opacity in the 
accretion disk, and m = 3/10, n = 4/5 when the free-free and free-bound transitions do. The "diffusion constant" D, 
defined by the specific vertical structure, relates Eq, F, and h: 



So 



2{l-m)Dh^- 



(7) 



(see also Filipov 1984 ). D is a function of a, opacity coefficient, and the dimensionless values, which are the combi- 
nations of the characteristic physical parameters of the disk. The value of D is to be derived from the consideration 
of the dis k vert ical structure. In the following section we obtain its value using the results of the work by Ketsaris & 
Shakura (|1998D . 



4. Vertical structure of standard disk 

Hereafter, until specially mentioned, we assume that the matter in the disk moves with the Keplerian angular velocity 
Wk, and its state is governed by the ideal gas equation 

P-^. (8) 

where ^ and 5R — 8.31 x 10^ erg mol^^K^^ are the molecular weight of the gas and the molar gas constant, respectively. 
Along the Z coordinate the hydrostatic equilibrium takes place: 

and the continuity equation is 

We assume the radiation transfer equation in the diffusive approximation: 
c 9(aT4) 



3 p dZ 



= -Q, (11) 



where c = 2.99 x 10^° cm s ^ is the light velocity, a — 7.56 x 10 erg cm'^K''. The vertical gradient of the radiation 
flux Q is proportional to the energy release per unit volume in the disk; that is, 

II =e [ergcm-3s-i] . (12) 

We take the opacities in the form of a power law k = K^p^ /T'^ where = v = 0, kq = 0.4 cm^g~^ if >Ct: ^ Kff and 
= 1, v = 7/2, }<o = 6.45 X 10^^ cm^K''/^g~^ if m:^ ^ Xt.. Generally, in the optically thick disks the energy release can 
be described as a power law of temperature and density (Tayler 1980). In a sense the calculation of the disk structure 
resembles the calculation of stellar internal structure. In the present study two cases are considered: the energy release 
£ is proportional to (a) the pressure (x pT, (b) the density p alone. The thermal energy release is due to the differential 
rotation of a viscous disk: 



a7 2' 



"np -jr- — 7:'^K Wrip ■ (13) 



We follow Shakura (|T97|) and Shakura & Sunyaev (|T97|) in suggesting that the turbulent viscous stress tensor is 
parameterized by the pressure: 

Wnp = - vt pr — = -uji^vt p = a P , (14) 
where vt is the kinematic coefficient of turbulent viscosity. The height-integrated viscous stress tensor is given by 

Wr^{r,t) = 2 Wr^dZ =3wk ivtpAZ. (15) 
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The energy emitted from the unit surface of one side of the disk is obtained by integrating (12) using and (|l|): 

Qo^l Wr^ [r, t) . (16) 

Above equations written for stationary accretion disks hold in a non-stationary case taking into account that 
the characteristic hydrostatic time of order of a/uj is shorter than the time of radial movement in the disk r/vr ^ 



There are now various works investigating the vertical structure of the disks. For example, Nakao & Kato (1995) 
considered turbulent diffusion in the disk providing the variations of viscous heating and a-parameter along the height 
Z . The vertical structure of the disks including radiative and convective energy transfer was investigated by Meyer 



& Meyer-Hofmeister (1982). They investigated two types of viscosity, proportional to the gas pressure or to the total 



pressure. We adopt the result of Ketsaris & Shakura (1998) who proposed a new method of calculating the vertical 



structure of optically thick a-disks assuming power p- and T- dependences for the opacity and the energy release. 
The dimensionless variable 

21] 

is introduced for convenient description of the problem, along with the following variables: -p = P/Pc, S = T/Tc, 
z = Z/Zo, j — p/pci and q = Q/Qo- The method involves the finding of the eigenvalues of the dimensionless 
parameters in the differential equations that describe vertical structure of the disk: |^ 

- n n ^ • n - ^° ^ • 
-^-nin^z, ni-^^. 



da p 2Zo Pc 

dq 3acjK3?TcSo _ a^Tcl^o 

da 4 QoP WripP 

n4 ■ U4, = — ( — ^ ^ ^° ^° 



(17) 



da 6"+^ ' 32 V / ' 

using the definite boundary conditions in the disk. Tc, pc, Pc denote the values in the equatorial plane of the disk and 
Qo = (ac/4)T4. ^ ^ ^ ^ . . 

After some algebraic manipulation of the right hand equations in (n^) we obtain So written in terms of Wr,p 
and u! r^, which in view of (|^) yields: 

D=^rr^\ „ (-1 (^] (GM)i^+«4 , (18) 



4(1 -m) \ Vm/ VSac 

where 

4 + 2C 12 + 11C-2I/ 

ITl — Tl 

10 + 3C + 2i^' 10 + 3C + 2z^ ■ 

It is worth noting that D depends on >ro very weakly: to a power of 1/5 or 1/10. This fact is believed to reduce the 
effect of uncertainties in our knowledge of the real law of opacity. The combination of ni_2,3,4 in ( |l8|) varies slightly 
with the optical depth r, i.e. along the radius of the disk (see Tables 0, §). Thus, factor D in the basic equation of 
time-dependent accretion ^ is considered to be constant. 

Specific energy dissipation e/ p — dQ/dT, is defined by the temperature variations over Z. In principle, the intensive 
stirring in the disk can account for the situation when the energy release per unit mass does not depend on the height 
Z. This refers to the case (b) mentioned above where e is the function of density. In this situation the temperature 
dependence disappears in the energy production equation (third line of (jlj)) and — 1. 



Ketsaris & Shakura (1998) calculated the values 01111,112,113, and 114. Selected values of ni.2,3,4 and corresponding 



7 /2 

values of S, in the Thomson opacity regime, and some effective optical thickness of the disk tq = So >co /9c/(2 Tc ), in 
the free-free regime, are shown in Tables |l|, |[ For the full version of Hi. 2, 3. 4 list and discussion the reader is referred to 

^ Left bottom equation in formula is corrected in comparison with the journal variant 
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Table 1. Vertical structure parameters in the Thomson opacity regime 
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Ha = 1 


log 5 


Hi 


Ha 




n4 


111 


ii2 


II4 


4.00 


6.37 


0.516 


1.150 


0.460 


6.46 


0.511 


0.500 


3.00 


5.67 


0.546 


1.149 


0.459 


5.74 


0.542 


0.499 


2.00 


4.47 


0.610 


1.142 


0.454 


4.52 


0.605 


0.490 


1.00 


2.61 


0.740 


1.105 


0.398 


2.63 


0.737 


0.417 



Table 2. Vertical structure parameters in the free-free opacity regime 



logro 


Hi 


Ha 


Ha 


Hi 


Ha = 1 


Hi 


Ha 


n4 


4.00 


7.07 


0.487 


1.131 


0.399 


7.12 


0.485 


0.437 


3.00 


6.31 


0.515 


1.131 


0.398 


6.34 


0.514 


0.436 


2.00 


4.98 


0.576 


1.126 


0.395 


4.98 


0.576 


0.431 


1.00 


2.83 


0.716 


1.095 


0.354 


2.81 


0.716 


0.373 



the original paper by Ketsaris & Shakura ( |1998D . The parameter S was introduced by them for the sake of convenience 
and denotes the ratio of total scattering optical thickness >Ct So to that at the thermalization depth|^; 



So/2 

Mr* = 1; 



{ng m:t^)^^'^ pdZ , 



(19) 



where r* is the effective optical depth (Zeldovich & Shakura 1969; Mihalas 1978). 



5. Time-dependent accretion in Keplerian disk 



5.1. Solutions to non- stationary Keplerian disk equation 

The self-similar solutions of Eq. (|^) were found by LS87. In these solutions any physical characteristic of the disk, for 
instance, the surface density Eo(f, i), can be presented in the form: So(r, i) = S{t) s{r / R{t)), where the scales S{t) 
and R{t) depend on t in a particular way, and s{r/R{t)) is a universal function of one self-similar variable ^ = r/R{t) 
(Zeldovich & Raizer 1967). The solutions represent three stages of the non-stationary accretion on an object. The first 
stage is the formation of the accretion disk from some finite torus around an object. The second stage is the developing 
of the quasi-stationary regime of accretion, and the third - the decay of accretion when the external boundary of the 
disk is spreading away to infinity. LS87 obtained the self-similar solutions of type II for the first two stages and the 
self-similar solution of type I (Zeldovich & Raizer 1967) for the final stage (when there is conservation of the total 
angular momentum of the disk). 

In a binary system the accretion picture has particular features. The main feature is the limitation of the outer 
radius. Thus, one cannot apply the LS87 solution at the third stage, that is, during the decay of accretion after the 
peak of the outburst. The spreading of the disk is to be confined by the tidal interactions. The tidal torque produced 



by a secondary star has strong radial dependence (Papaloizou & Pringle 1977). As Ichikawa & Osaki (1994) showed, 
the tidal effects are generally small in the accretion disk, except near to the tidal truncation radius, which is given 



by the last non-intersecting periodic particle orbit in the disk (Paczyhski 1977). They concluded that once the disk 
expands to the tidal truncation radius, the tidal torques prevent the disk from expanding beyond the tidal radius. 

A class of solutions of Eq. (^, which this paper focuses on, can be found on separating the variables h and t. We 
seek the solution in the form F{h,t) — F(t)/(^), where ^ — h/ho, ho — (GMrout)^^^ ■ From Eq. (|^), substituting 
h = h^, it follows that 



Mih,t) = ~2TTf'{h/ho)F{t)/K ■ 

Substitution of the product of two functions in Eq. 

un+2 \ 1/"' 

F{t) ■ 



gives the time-dependent part of the solution: 



-XmD{t + to) 



(20) 



(21) 



Formula (|l|) is corrected in comparison with the journal variant 
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Fig. 1. The solution /(^) in two cases: when x.-^ ^ (solid line) and >fff 3> >Ct (dashed line) 
Table 3. Summary of parameters in solutions for two opacity regimes for the Keplerian disk 





m 


n 


A 


ao 


ai a2 


k 




yii > xft 


2/5 
3/10 


6/5 
4/5 


-3.482 
-3.137 


1.376 
1.430 


-0.39 0.02 
-0.46 0.03 


3.8 
3.5 


6.6 
6.0 



D is the constant defined by the vertical structure of the disk (§ 4, Eq. (|18D); A is a negative separation constant which 
can be found from boundary conditions on /(^); to is an integration constant. From here on we set io = for the 
Thomson opacity regime. We calculate a value of to for the free-free opacity regime in § 6.2. Expression ( pT| ) represents 
asymptotic law for after-peak evolution of a real source. 

The equation for /(^) is a non-linear differential equation of second order which is a particular case of the general 
Emden- Fowler equation (Zaycev & Polyanin 1996): 



" J _ MU rl- 



(22) 



the solution of which we seek as a polynomial 

/(O - "0^ + + 02^' + . . . . 

Substituting /(^) into Eq. (02) we obtain for the second and the third term: 



fc = 3 



/ = 2fc- 1, 



TO, ai 



a2 



k{k - 1) 



(23) 



(24) 



1(1-1) 



(1 



and ao, A are to be defined from the boundary conditions on /(^). 

We consider the size of the disk to be maximum and invariant over the period of outburst. As the drain of angular 



momentum occurs in a narrow region near this truncation radius (Ichikawa & Osaki 1994), we treat the region near 
this radius as the (5-type channel, not considering the details of the process. In other words, the smooth behaviour 
of spatial factor / in the moment of viscous forces F (which increases as oc r^/^ in the inner parts of the disk, then 
flattens, reaches the maximum and drops down near rout due to tidal torque) is analytically treated as increasing, 
flattening, and reaching maximum at rout, which is the end of the disk (this profile is shown in Fig. |l|). Thus we 
propose the boundary conditions as follows: 



/(1) = 1, /'(1) = 0. 



(25) 



Corresponding ao and A are displayed in Table |[ 

Naturally, real accretion disks have finite value of rjn ^ 0, but still, in most cases, rin/rout ^ 1, that is equivalent 
to rin/rout = in our problem from the mathematical standpoint. 

Note that ( pi] ) implies a considerably steeper time dependence than the solution by LS87 does. The latter yields 
the accretion rate as a function of t-i^/i^ jf -^^ ^ -^^^ a,nd if >cg ^ xtt- In our case this dependence is t~^^^ if 

>Ct. ^ >Cf{, and t~^^/^ if 3> >Ct- This difference is due to the non-conservation of angular momentum in the disk in 
our case. 

The following subsections contain the explicit expressions for the physical characteristics of the disk. They are 
deduced from (Q), (p7|), (|l8|), and (|2l|). We introduce for the mass of the central object the quantity = M/Mq. 
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5.2. Thomson opacity regime (>Ct- >Cff) 

Here the function / = /(^) = /((r/rout)^^^) and the values Hi. .4 should be taken for the Thomson opacity regime. 
Then we have^: 



^ [g-2/5^j^28/5g-i7/5] = 2.42 X lO^^ a^/^ m^/^ (-^) '^'(n4n4 
So [gcm^^] = 3.4 X 10^ a~^m^^'^ ( I 



0.5/ V^out 

11/10 



-9/10 



/ 



3/5 



/ \ 3/2 / 



/ \ —i-i-/i-^> / \ 1/2 / 

T.|K|.l,8xlO.„-.„;/.(ii)(j_) ( 

/ \ / + \ -1/2 



-3/2 



1/2 /^x -1 



-1/20 



/2 



T* = 4.8 X 10^ a-^ <' ^ 



5/4 / r \i/io 



.1/10 (^rsiiiV^' f^y'^' 

0.5/ Vwy v^©; Viov Vninl 



/4 



(26) 
(27) 
(28) 
(29) 
(30) 



r* is the effective optical thickness of the disk defined by the combined processes of scattering and absorption. We 
take approximately (c.f. (|l9|)): 

^ _ ^^ 0.4 X 6.45 X 10^^ pc ^ 



7/2 



5.3. Free-free opacity regime (>Cf[ 3> x^tJ 

Here the function / and the values Hi. .4 should be taken for the free-free opacity regime. The following formulae 
contain the constant to appeared in expression (^l|). It was neglected in the previous subsection; here to accounts for 
the possibility of time shifts between the solutions in the two opacity regimes. We have^: 



D [g-3/io s-16/5] ^ 5.04 x lO^^ a^/s (^ii^ ^8 ^ 



-1/10 



Eo [gcm-2] = 5.3 X lO^a-s/^mye (J^y^^ ( J_ 

V0.5/ \rout 



-11/10 / X 13/6 

^7/10 



\Rq) \ lOd 



-10/3 



Tc [K] = 3.1 X 10^ a" 



/ X -9/10 / \ 1/2 / , , , \ 



-1/2 



(HiHa)^/' , 



(31) 

(h}/'H2H|H4)'^'(32) 
(33) 
(34) 



r = 4.5xl0^a-V3^i/6 ( ' 



\ront 



HI Hi 



n}/^H2, 



(35) 



This regime is characterized by lower temperature and density, and the optical thickness of the disk is defined by the 
processes of free-free absorption: r = 2to = 6.45 x 10^^ pc Tc "^^^ Eo . 



Formulae (2£) and 



29) are corrected in comparison with the journal variant 



Formulae (31) and (34) are corrected in comparison with the journal variant 
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6. Bolometric light curves of time-dependent standard accretion disk: power law 

In order to calculate the luminosity of the disk, we assume the quasi-stationary accretion rate as it is at r <C ^out- 
For these most luminous parts of the disk we take M(t) = M{0,t) given by (pO|). The overall emission of the disk is 
defined by the gravitational energy release L — rj M{t) (? , where rj is the efficiency of the process. 

At early t, when the Thomson scattering is dominant, we derive that the bolometric luminosity of the disk varies 
as follows: 

= 1.3 X 10» o-„4" (^) {-^S {^f' {4^y"' (mn,)"' . ,36) 

As the temperature decreases, the law of decline switches to: 

L,it) [ergs-]=3.6xl0-a-/3^^/^ (il) 0^)^^'^ {^f'^ (^)"°'' (n}/^ nl H,) . (37) 
The mass of the disk can be derived by integrating So over its surface: 

Mdi.k,T [M0] = 4x lO-V-^mV^ (iL)' {"^^ (l^) i^t^^f^ (38) 

Af...,. [AfQ] = 2 X io-«-/3™V6 ^±_y/' (^)"'" (^)"'' (n;/^n.n«n,)^^^ . (39) 

The constant tp is the same as in the previous section. These solutions give an asymptotic law for the disk bolometric 
luminosity variations. The value of to will be obtained in § 6.2 when we shall discuss the transition between the regimes 
of opacity. 

We remark that the observed X-ray light curves can have different (most probably, steeper) law of decay. Indeed, 
the energy band of an X-ray detector usually covers the region harder 1 keV where the multi-color photon spectrum 
of the disk (having appropriate temperature) can have turnover from —2/3 power law into exponential fall. This 
turnover is expected to change its position due to variations in temperature of the disk after the burst. The narrower 
the observed band, the more different the observed curve could look like in comparison with the expected bolometric 
flux light curve. In § 7 we discuss this subject in more detail. 



6.1. Luminosity - accretion disk parameters dependence 



It is essential to point out that in formulae (|36D, ( p7D the parameters a, fi, m, rout, Hi, 2, 3, 4 cannot be changed to describe 
how luminosity depends on them. Indeed, these expressions were found as a result of solution of differential Eq. ^ 
with the constant coefficient D (which depends on parameters of the disk, except 77). Imagine a situation when one of 
these parameters, say a, quickly increases. This will not result in the decrease of the luminosity as it might seem from 
( p^ ) or (^7|). What will happen really is that the accretion will change to another solution (during the same regime of 
opacity), according to the new D*. Supposing that the mass of the disk remains constant during this transition and 
taking into account that the profile of So(r) does not change, it can be seen from (^ that F changes discontinuously 
and the luminosity L oc M oc F{t) jumps as (a*/a)''/'^'^^~™'-', a* > a. The relation between the new and the old F 
and D, obtained from (|7|), gives the new term (t + tQ) in (pl|): 



t + t* 



t + to \D 



(40) 



Thus the increase in a gives the increase in D and, consequently, (i-t-tg) < (* + ^o) which implies a steeper light curve 
after the transition than before. The increase of a can be possibly provided by the enhanced role of convection in 
the accretion disk and will result in the brightening of the disk. This situation is displayed in the inset in Fig. |^. We 
note that the descending portion of the curve after the increase is uncertain if convection is involved since the disk 
structure modifies from that presented in § 4. 
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6. 2. Thomson opacity - free-free opacity transition 

The temperature of the disk decreases with time, and eventually the free-free and free-bound opacity supersedes 
the Thomson one. It is possible to connect two regimes at the point ttr), where Fi(^,ftr) = ^2(^1 ^tr + ^o) and 
Eo.i — So_2 (indexes 1, 2 denote different opacity regimes) - two conditions allowing us naturally to define both ttr 
and tQ-. 

/ // \^/^ / r \ -4/5 /„ \4/5 fri2/5 /fn fn2 fni6 fn2 \ 

(41) 



The right top indexes of / = /(^) and Hi 2,3,4 indicate the opacity regimes. As the profiles of /(^) = /((f/rout)^^^) are 
very close in these two regimes (see Fig. and parameters 1X1,2,3,4 vary slightly with radius (being roughly constant 
in the region where the substantial mass of the disk is enclosed), the physical parameters of the disk {'Soir), Tc(r), 
etc.) calculated in the two solutions are sufficiently accurately equal. 

At the time ttr the free- free absorption coefficient xs — xop'^/T^ calculated in the Thomson opacity regime and 
in the free-free opacity regime takes the form: 



f r^„,2„-ii _ n Qoo 



'>fff [cni^g-^] = 0.399 



-..ff [cmV^] = 0.399 (^^ j 



'f\' fn|fn4 



(42) 



The closeness of >cg to = 0.4 cm^g^^ confirms the reliability of our calculations and yields the smoothness of the 
transition. 

Fig. H represents the bolometric light curve of the disk for a — 0.3, = 3. Hereafter we substitute ni_2,3,4 with 
their typical values in a self-consistent way. The transfer between Thomson and free-free regimes begins at the moment 
r/rout = C" = 1, < « 8'^(tox/3)2/^ a'^/s {^/0.5f/^ (rout/i?©)^/^ - arrow A at 21^^ in Fig. |. We intersect the curves at 

r/rout = 0.5, t = ttr ~ 13'^(m,/3)2/5 a'^/' i^i/O.bf/' (rout/i?©)"/' , 

what corresponds to itr ~ 34'^ and to ~ 17'^ for a = 0.3 (left small arrow). We call ttr "moment of transition". The 
transition ends at the time t « 21'^(mx/3)^/^ a^^/^ {rout/ Rq)^^^ when the solutions match at r = 0.25 rout - 

arrow B at 55*^ in Fig. ^j. 

This picture is reliable and useful, even though it implies the existence of two separate regimes, which is evidently 
not quite true. Indeed, at any epoch the inner part of the disk would be scattering dominated, the lower the accretion 
rate, the smaller this part. Obtaining of an exact solution needs consideration of combined free-free and Thomson 
opacity of the gas. 

There is some t which corresponds to the Eddington limit L-^^d ~ 1-3 x 10"^^ erg s^^. This means that the 
real source evolution could be described in our model only at later t. Thus, generally speaking, the solution before 
this moment appears inapplicable. As seen in Fig. ||, the applicable part of the solution belongs almost entirely to the 
free-free opacity regime (the bold dashed line). 

The second intersection of the curves in Fig. || at t « 95'^ (right small arrow) corresponds to the other intersection 
of functions Fi(^,t) = F2{£,,t + to), meanwhile the physical parameters of the disk calculated using formulae (^6|)--(^5|) 
are different. Thus the disk is at the same (free-free) opacity regime as before. 

When Tc decreases to the value ~ 10^ K, the convection (which presumably appeares in the zones of partial ioniza- 
tion) starts to influence the disk's structure, and the diffusive type of radiation transfer, which we use, is no longer valid. 

For TOx = 3 and a = 0.3 this happens at t « 190^' : t+to « 2,2'^ rnj^ a-^ (Ai/0-5) (rout/i?©W" (r/rput)''^/ ^" ^ f^ '^° ^^s. ■ 
For inves tigation of the disk evolution on larger time-scales see e.g. Cannizzo et al. ( 1995|) , Cannizzo ( 1998 ), Kim et 
al. ( 119991) . 



7. Observed light curves 

As we mentioned in § 6, the observed light curves can have a slope of decline which is different from that of the 
bolometric light curves due to particular spectral distribution. In this section we are going to illustrate this suggestion 
assuming the simplest spectral distribution of the disk emission. 
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50 100 150 200 

t, days 

Fig. 2. Bolometric luminosity Lt and Lg calculated for parameters: = 3, a = 0.3, = 0.5, Tout = ^0- Shown 
are the solution in Thomson opacity regime (solid line) and in the free- free opacity regime (dashed line). Their bold 
parts represent the resulting light curve of the disk. Small arrows mark two intersections when Fi (^, t) = F2{£,,t + to). 
The inset illustrates the case of increase of a from 0.3 to 1 




To calculate the spectra, one can assume the quasi-stationary accretion rate in the inner parts of the disk because the 
M variation is small there (M oc f '{£,), see Fig. |l|). The outer parts of the disk, where accretion rate varies significantly, 
contributes to the low-frequency band of the spectrum. We discuss the X-ray band and the most luminous parts of 
the disk and, thus, we take M{t) = M{0, t) given by (f^). 

Provided each ring in the disk emits as a black body, the temperature of the ring can be found as follows: 



1 d 3 

<J5BT^^^-Wr^r^ = - C^K Wr^ , (43) 
where ctsb — 5.67 x 10^^ erg cm^^s^^K^^ is the Stephan-Boltzmann constant. Then 




1/4 

y Sttctsb^"^ L \ r ) J 

In the last expression the stationary solution for Wnp is taken. The black-body approximation is satisfactory if 3> >fT- 
Then the outgoing spectrum is the sum of Planckian contributions of each ring of the disk and has the characteristic 1/3 
slope for photon energies <C fcTmax, where Tmax is the maximum effective temperature of the disk (Lynden-Bell 19691 ). 



However, if the Thomson scattering on free electrons contributes substantially to the opacity, the outgoing spectrum 



is modified (Shakura & Sunyaev 1973). See e.g. Ross & Fabian (1996) for investigation of spectral forms of accretion 
disks in low-mass X-ray binaries. 

The light curve is simulated by integrating at each t the spectral density 

rout 

^ K f rdr 

c2 J exp(/ipiy/fcr(r,t))-l ^ ' 



over the specific frequency range using ( pO| ) and (Q), where hp = 6.626 x 10~^^ erg s is the Planck constant. The 
numerical factor in ( p5[ ) corresponds to the luminosity outgoing from one side of the disk. 

Explaining the observed faster-than-power decay of outbursts in soft X-ray transients, one must take into account 
the specificity of the energetic band of the detector. Naturally, the observed slope of the curve depends on width and 
location of the observing interval. The narrower this band, the more different the observed curve could look like in 
comparison with the expected bolometric light curve. Of course, this difference also reflects the spectral distribution 
of energy coming from the source. 

We show here how the slope of the curve changes in the simplest case of multi-color black body disk spectrum 
according to which spectral range is observed. Following ( ^5|) we calculate and integrate it over three energy ranges: 
3-6 keV, 1-20 keV and that one in which practically all energy is emitted. Fig. || shows the photon flux variations in 
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Fig. 3. The flux from one side of accretion disk at 1 kpc for parameters: mx — 3, a — 0.3, fi — 0.5, Tout — Rq- The 
curves show the bolometric flux (upper curve), the 1-20 keV flux (middle curve) and the 3-6 keV flux (lower curve) 
during the Thomson opacity regime (solid parts) and the free-free opacity regime (dashed parts) 




two X-ray energy ranges (those of Ariel 5 and EXOSAT oi Ginga observatories) and the bolometric flux variation for 
the face-on disk at an arbitrary distance of 1 kpc. The vertical line marks the time after which bolometric luminosity 
of the disk's one side is less than LEdd- 

One can see an almost linear trend of the X-ray flux when bolometric luminosity is under the Eddington limit (to 
the right of the vertical line in Fig. especially in intervals of ~ 50"^. The decline becomes closer to the exponential 
one with time. The slope of the curve depends on a, m, rout) and other parameters. For the same parameters as in 
Fig. ^, the e-folding time falls in the range 20-30 days for the lower curve (3-6 keV). For instance, smaller a will result 
in less steep decline. 

The natural explanation of such a result is the following: because the spectral shape of the disk emission has Wien- 
form (exponential fall-off) at the considered X-ray ranges, the law of variation of X-ray flux is roughly proportional 
to exp^-hpf/kT^^it)). In the free-free regime of opacity we have T^^it) cx ^^^''(i) oc M{ty/* oc ^-10/12 (gee § 5.1). 
Consequently, the observed X-ray flux varies like exp (— t^/^), which is quite close to exponential behavior. We restrict 
ourselves to this brief and general discourse, as a detailed application of our model to observed sources is not a goal 
of this paper. 



8. Viscous evolution of advective disk 



As we know, the structure of an accretion disk in the vertical direction, the relation between the viscous tensor and the 
surface density in particular, defines the type of its temporal evolution. In advective disks, which are the low-radiative 
accretion flows, the relations between their characteristic physical parameters differ significantly from those in standard 
disks. In this secti on, we discuss the re sults of § 2 as applied to the disks which radial structure was presented by 
Spruit et al. ( |l987D and Narayan & Yi ( |1994 |l995i hereafter NY). 

The viscous stress and the surface density are related through the kinematic coefficient of turbulent viscosity. 
Integrating the component of viscous stress tensor one obtains (c.f. ( p^ and (p^))p|: 



Wr^{r, t) = 2 I Wrcp dZ = -2 j pvt^ r dZ = - ^ rY^^Vt , 



(46) 



where vt is the averaged kinematic coefficient of turbulent viscosity. Then the relation between T,o{h, t) and F = Wnp 
is given by 



dh 



(47) 



Formula is corrected in comparison with the journal variant 



12 



G.V. Lipunova & N.I. Shakura: New solution to viscous evolution of accretion disks in binary systems 



Recall that /i* is the real specific angular momentum and h is the Keplerian one. It can be seen that i>t{h, t) and 
define what class of solutions Eq. (^) will have. 

If one adopts for the structure of advection-dominated accretion flow (ADAF) the self-similar solution by NY, it 
can be easily inferred that such disks exhibit the exponential with time behaviour. The solution of NY is given by: 

(48) 



Vr = -Ci Wk r , OJ ^ C2UJk , 

Expressing Eq in the basic Eq. 
dF _i _ {GMf d^F 

Solution (E 



2 2 

C3 w„ r 



in terms of F, we obtain from M7t) and (48): 



enables deriving the relation between Pt and h 

du) 3 _ 2 
-g;^ ^ I't P C2 uJk = a p , 

where Og is the isothermal sound speed. Thus i>t is a function of radius alone 
vt 2 a flg 1 



lok r"^ using the a prescription of viscosity: 



vtpr- 



2 C3 

- —a 

O C2 



3c2Wk WkT^ 

and Eq. ( |46| ) can be rewritten in the form: 
dF _ Da d^F 
'dt ^ 



(49) 



(50) 



(51) 



h dh^ 



{GMf 

2C2 



(52) 



h/ha, ho being some value of h: 



.2 =^e./ 



(53) 
(54) 



Solution to (52) is sought as a product of two functions /(^) and F{t), with ^ 
F{t)^F° e:>^p{XD,t/hl) , 

The exponential temporal behaviour of NY flow is evident. Generally speaking, any disk possessing such properties of 
vt as constancy in time would have such exponential behaviour because its evolution would be described by a linear 
equation (like (p2[)). 

The question is, would the confined NY disk keep such properties or it would not. The fact is that NY solution 
describes the infinite disk. Either the boundary conditions destroy the linearity of (p3) or just the characteristic decay 



time changes, this problem requires further accurate numerical investigation. For instance, Narayan et al. (1997) cal- 
culated numerically the global structure of stationary advection-dominated flow with consistent boundary conditions; 
they noted that although the self-similar solution (48) makes significant errors close to the boundaries, it gives the 
reasonable description of the overall properties of the flow. 

Further we assume that exponential trend of solution persists. Generally speaking, the equation determining f{£) 
will differ from (|4|). This difference may be not very significant. One can see that E g. ( I s^) is a particular case of (E2|) 
where n — \ and m — Q and, hence, the solution can be found according to (|2^ ) and (^). Besides, the solution of (54) 
can be found in terms of Airy functions (Bcsscl functions of order 1/3). 

The accretion rate evolves with time as follows (c.f. ( ^o|) and (|5^)): 

Ai^ ocexp(AL>a</^o) • (55) 

The value of accretion rate can be determined if an initial condition is imposed at some t. Mahadevan ( 1997 ) showed 
that ADAF luminosity oc Af^ or ex; M according to whether the electron heating is dominated by the Coulomb 
interactions or by the viscous friction. Subsequently, the luminosity has an exponential decay too. 

We can estimate the characteristic time of evolution of such flow. It can be obtained from (pq). Let us compare 



the diffusion time 

tov 1 Cn 



r^/z/t with the corresponding orbital period 27r/u;. Since vt 



aag/w, with ([48|) we have: 



1 



5/3-7 



torh a 2 7rc3 2 7ra(7-l)/a ^^^^ 
We use the expressions for C2 and ca from NY. Here 7 is the ratio of specific heats; /a measures the efficiency of 
radiative cooling. In the limit of no radiative cooling, we have /a = 1 while in the opposite limit of very efficient 
cooling /a = 0. NY solution is degenerate if 7 = 5/3 because the angular velocity of the flow is zero in this case. 

We can see that the time-dependent advection-dominated disk is quickly depleted if a is not small. For example, 
consider the light curves of X-ray novae which have the exponential decay time scales ~ 30"^. To obtain tev of such 
order, a should be ^ 10^^. However, the advection-do minat ed solution ceases to exist if the accretion rate is greater 



Edd 



than the critical value a 

1.39 X 10"'^^ m^. Hence, a ~ 10^^ yields the critical accretion rate M, 



(Narayan & Yi |1995| , Mahadevan |1997| ), where MEdd = ^Edd iv/O-l)' 



10' 



g s" 



Edd- 
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9. Discussion and conclusion 



In this work, we presented the analytical solutions to time-dependent accretion in binary systems. During an outburst 
we propose the specific external boundary conditions on a disk confined due to tidal interactions. For two opacity 
regimes the full analytical time-dependent solutions for the Keplerian disk are obtained and an asymptotic light curve 
is calculated with smooth transition between opacity regimes. During the decline phase accretion disks around black 
holes appear to be dominated by the free-free and free-bound opacity in order to comply with the Eddington limit on 
luminosity. This phase is characterized by the power-law decay of accretion rate oc i^^"/-^. It is shown that the decay 
time scale depends on the real energetic band of detector (Fig. ^ . 

The results obtained in this work can be applied to the accreting systems having variable emission of flare type if 
emission is essentially due to the fully ionized accretion disk around a black hole, or a neutron star, or a white dwarf. 

Narayan & Yi (1994) accretion flows are shown to undergo exponential decays if the disk has infinite size. This 
notable result probably persists even when the advective disk is in a binary system. The latter suggestion is to be 
thoroughly considered in the accurate numerical investigation. If this is the case, the abrupt steep falls observed in 



several X-ray novae (Tanaka & Shibazaki 1996) in the last phase of the decay, at luminosity levels < 10 erg s , can 



be interpreted in terms of quickly depleting ADAF (§ 8) with relevant values of a ~ 10^ . 

Using the results of this work, we can explain the general features of X-ray novae light curves in the early phase. 
Typical XN outburst light curves (see Tanaka & Shibazaki [1996 Chen at al. 1997 for a review) show quasi-exponential 
decay. To date several approaches have been used to accoun t for X N featur es. Th e exponen tial d ecays were obtained 
in the framework of disk instability model (Cannizzo et al. 1995; Vi shniac 1997 ; Cannizzo 1998 ) in which the large 
time-scale evolution of the disk is considered. Mineshige et al. ( 1993 ) argued that the exponential decays in XN can 
be reproduced if the mass and the angular momentum are efficiently removed from the inner portio ns of t he disk at 
a constant rate, or wind mass loss or enhanced tidal dissipation could be substantial. King & Ritter (1998) took into 
account the irradiation of the disk by the central X-ray source and obtained the characteristic XN light curves. 

We suggest an alternative reason to explain this remarkable feature, at least during the early stages of the outburst 
when the disk is fully ionized. Nearly exponential X-ray decays oc exp {—t^^^) are obtained taking into account the 
fact that the X-ray light curves are observed in the energetic range where the spectrum of the disk has Wien-form. 
Black hole XN spectra typically are composed of an ultrasoft component and a hard power-law component (e.g. 
Tanaka 1992; Tanaka & Shibazaki 1996). At the first stages after outburst the ultrasoft component dominates and can 



be represented by a multicolor blackbody disk (Tanaka 1992). This component has an exponential fall-off, a decisive 
factor to produce observed exponential trends. The observed characteristic times can be obtained within reasonable 
intervals of parameters (Fig. ||, § 7) . 

The secondary peak commonly observed in XN can be qualitatively analytically produced by certain reconstruction 
of viscosity mechanisms and corresponding increase of a (§ 6.1). Possible mechanisms of reflares involving irradiation 
effects were investigated by Kim et al. (1994), Mineshige (1994), King & Ritter (1998) (see, however, Cannizzo 1998). 

Of course, the accretion disk spectrum represents only one contribution to the total observed spectrum of the 
source. The corona around the disk is probably responsible for the other spectral components. In addition, taking 
into consideration the irradiation of the outer parts of the disk would affect evolution of the disk (see, e.g. King & 

1999). Kim et al. ( |1999D constructed an optical light curve of a XN and found the direct 



Ritter 1998: Kim et al 



irradiation of the disk by the inner layers to have only a small effect on the outer disk because of shadowing. The 
indirect irradiation (from a corona or a chromosphere above the disk) is found to affect the light curve more strongly. 
We suggest that the irradiation of the twisted warped disk could also result in important heating of the outer layers. 
Further investigation and applications to observed sources will be the basis of our future work. 
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